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Abstract We study dynamic pricing mechanisms of European contingent claims under 
uncertainty by using G framework introduced by Peng (|21 1). We consider a financial market 
consists of a riskless asset and a risky stock with price process modelled by a geometric 
generalized G-Brownian motion, which features the drift uncertainty and volatility uncertainty 
of the stock price process. A time consistent G-expectation is defined by the viscosity solution 
of the G-heat equation. Using the time consistent G-expectation we define the G dynamic 
pricing mechanism for the claim. We prove that G dynamic pricing mechanism is the bid-ask 
Markovian dynamic pricing mechanism. The full nonlinear PDE is derived to describe the bid 
(resp. ask) price dynamic of the claims. Monotone characteristic finite difference schemes for 
the nonlinear PDE are given, and the simulations of the bid (resp. ask) prices of contingent 
claims with uncertainty are implemented. 

Keywords G Brownian motion, G expectation, uncertainty model, dynamic pricing mecha- 
nism, monotone finite difference 



1 Introduction 

In probability framework, the uncertain model of the stock price is assume that the stock price be a 
positive stochastic process that satisfies the generalized geometric Brownian motion 

dlogS t =ijfdt + ofdW t , (1.1) 
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where (W t ) t >o be a 1 -dimensional standard Brownian motion defined on a probability space (£2, JF, (!ft),P) 
and (^F, : < f < P) is the filtration generated by W ( , and (/jf ,af ) r >o is unknown such that 

(/jf ,of ) G [/J,ju] x [a, a], /j,/7,o,a are constants such that /j < /7,o < a. (1.2) 



We denote T as all possible paths (/jf , of ),>o satisfying ( 1 .2 ), then T is a closed convex set. For each 
fixed path (jt)t>o = (j^>^t)t>0 G T, let Py be the probability measure on the space of continuous 
paths (C(0,oo),(8(C(0,oo))) induced by Jd^f^ + of t/W s ), and denote Po as the reference probability 
measure induced by W t . We set T be the class of all such probability measures Py, and for each Pef 
we denote the corresponding expectation. 

If the uncertainty comes from fjf which is called drift uncertainty, each probability P G fP is 
absolutely continuous with the reference measure Po. Chen and Epstein [7| proposed to use g- 
expectation introduced in [18] for a robust valuation of stochastic utility under drift uncertainty. 
Karoui, Peng and Quenez in ifTSIl and Peng in |[T9l proposed to use time consistent condition g- 
expectation Ef T : L 2 (j7r) — > L 2 (y t ) defined by solutions of a BSDE, as bid-ask dynamic pricing 
mechanism for European contingent claims £, 

ask price: Ef T [H' T Q = sup£p[5^| %\, bid price: -Ef^-H 1 ^] = - sup Pe2 ,£p[-H^|J f ], 
/■ ■/■ 

where (Hl;s > t) be some deflator. Delbaen, Peng and Rosazza (|8|) proved that any coherent and 
time consistent risk measure absolutely continuous with respect to the reference probability can be 
approximated by a g-expectation. 

The volatility uncertainty model was initially studied by Avellaneda, Levy and Paras (2] and 
Lyons [ 14 1 in the risk neutral probability measure, they intuitively give the bid-ask price of a Euro- 
pean contingent claim as follows 

askprice: sup E Q [e-< T -%\%\, bidprice: - sup eG o E Q [e- r ( T -%\F t }. (1.3) 
2GQ. 

where Q, = {Q : Q is the risk neutral probability measure of P, VP G T}, r is the short interest rate. 
For the volatility uncertainty, the probabilities in T are mutually singular, the upper-expectation can 
not be approximated by any g-expectation. 

Motivated by the problem of coherent risk measures under the volatility uncertainty ([1]), Peng 
( ||2TI . ||221 ) introduced a sublinear expectation on a well defined sublinear expectation space, under 
which the canonical process (B t )t>o is defined as G-Brownian. The increments of the G-Brownian 
motion are zero-mean, independent and stationary and N({0},[g 2 t,a 2 t]), and the corresponding 
sublinear expectation is called G-expectation. By using quasi-sure stochastic analysis, Denis, Hu 
and Peng in [ 1 1 1 and Soner, Touzi and Zhang in E4ll constructed consistent G-expectation, condition 
G-expectation and G-Brownian motion such that the condition G-expectation can be represented as 
an upper condition expectation 

£, G H=esssup eG(i £ e g|^], lje4, Q-a.s. (1.4) 

and M t is a G-martingale if and only if 

M s = essswpQ e qEQ[M t \!F s ], 0<s<t, Q-a.s. (1.5) 

which means that the condition G-expectation could be used as bid-ask dynamic pricing mechanism 
under volatility uncertainty, and a G-martingale is a supermartingale for each risk neutral probability 
measure space. 

There is a literature on the pricing the contingent claims under volatility uncertainty recently, 
Denis and Martin lfl2l studied the superhedging of the claims by using quasi-sure anyalysis, and 
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Vorbink |27| derived the lower arbitrage price and upper arbitrage price based on the approach of 
arbitrage by using G-framwork. For more general situation, if the uncertainty comes from drift 
and volatility coefficients, Peng |20| studied the super evaluation of the contingent claims by using 
filtration consistent nonlinear expectations theory. 

In this paper we study the dynamic pricing mechanisms of European contingent claims written 
on a risky asset under uncertainty by using G-framework introduced by Peng (2005, [21 1). At 
first, in a path sublinear space (Q.,Lj p (Q.),E) we model the price process of the risky asset by a 
generalized G-Brownian which describes the drift uncertainty and the volatility uncertainty of the 
stock. Using the techniques in G-framework we show that the risk premium of the stock is maximum 
distributed. We define the bid-ask dynamic price of the claim by using BSDE and derive a bid- 
ask dynamic pricing formula by using G-martingale representation theorem l25l . Further more we 
define a time consistent G-expectation E G by the G-heat equation and define a corresponding G- 
Brownian. By the G-expectation E G transform, the uncertainty model is transfered to volatility 
uncertainty model in (Q.,Li p (p.),E ), we prove that the condition G-expectation E G T is the bid- 
ask dynamic pricing mechanism for the claims, we also show that the bid-ask dynamic pricing 
mechanism E' t T is a Markovian dynamic consistent pricing mechanism and characterize the bid 
(resp. ask) price by the viscosity solution of a full nonlinear PDE which is the Black-Scholes- 
Barenblatt equation intuitively given by Avellaneda, Levy and Paras([2|). For numerical computing 
the full nonlinear PDE, we propose monotone characteristic finite difference schemes for discrete 
solving the nonlinear PDE equations, provide iterative scheme for the discrete nonlinear system 
derived from the characteristic difference discretization, and analysis the convergence of the iterative 
solution to the viscosity solution of the nonlinear PDE. In the end, we give simulation examples for 
the ask and bid prices of contingent claims under uncertainty. 

This paper is organized as follows. In Section 2, we give the financial market model. Section 
3 we derive a bid-ask price formula for the European contingent claim. In Section 4, we propose a 
G-martingale transform and G dynamic pricing mechanisms for the claims. Section 5 we investigate 
the Markovian case for the G dynamic pricing mechanisms and full nonlinear PDEs are derived to 
describe the bid and ask prices for the claims. Numerical methods for the nonlinear PDEs are given 
in Section 6. The implementations of the simulations for the digital option and butterfly option under 
uncertainty are shown in Section 7. 

2 The market model 

We denote by £2 — C(R + ) the space of all R— valued continuous paths (co r ), GR + with coo = 0, 
equipped with the distance 

p(co\co 2 ) := £2- i [(max Ico, 1 - cof | ) A 1], (2.1) 

then (£2, p) is a complete separable metric space. For each fixed T G [0,°°), we denote Q.J = {co.Ar : 
co G £2}. 

For the canonical process (Z? r )(co) = co ( ,f € [0,°°), for CO <G £2, we set 

L ip {Q. T ) {ty(B ti , • • • A,) : Vn G N,h,--- ,t n G [0, T], V<|> G C b , Up {R n )}, 

and 

L ip (£l) := U" =1 L ip (£l n ), 

where Cb±jp{R") denotes the linear space of functions cp satisfying 

|<PW -<P(y)| < C(l + \x\ m + \yD\x-y\ forx,y G R, 
some C > Q,m £ N depending on cp. 
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Assume that p., p, o and o are nonnegative constants such that p < p and o < o, we denote 
(£l,Lj P (£i),E) as a sublinear expectation space such that the canonical process (B t ) t >o is a general- 
ized G-Brownian motion with 

-E[-B t ]=pJ, E[B t ]=pt, 
-E[-B 2 ]=g 2 t 1 E[B 2 ]=o 2 t, (2.2) 

and the generalized G-Brownian motion can be express as follows 

B,=B,+bt, 

where (B t )t>Q is a G-Brownian motion andB, is A^({0}, [o 2 t,a 2 t]) distributed, and/?, is NQ/Jt ,pt], {0}) 
distributed. (Peng in [23 1 gave the construction of the sublinear expectation space (£l,Li p (£i),E), the 
generalized G-Brownian motion (B,) t >o, G-Brownian motion (B t )t>o an< iN(\pt,pt],{0}) distributed 
(b,) t > .) 

In this paper, we consider a financial market with a nonrisky asset (bond) and a risky asset (stock) 
continuously trading in market. The price P(t) of the bond is given by 

dP(t)=rP{t)dt, P(0) = 1, (2.3) 

where r is the short interest rate, we assume a constant nonnegative short interest rate. The stock 
price process S t solves the following SDE 

dS, = S,dB, , ( or dS, = S, (db, +dB,),) (2.4) 



where B t is the generalized G-Brownian motion. The properties showed in (2.2 1 imply that the 
generalized G-Brownian describe the drift uncertainty and the volatility uncertainty of the stock 
price. 

Remark 1 The generalized G-Brownian motion can be characterized by the following nonlinear 
PDE: u(t,x) = E[(p(x+B t )] is the viscosity solution of 

d t u(t 1 x)-g(d x u(t 1 x))-G(d xx u(t,x)) = 0, u\ t =o = q>(x) 

where <p(x) is a Lipschitz function, g(a) = ]JGL + — /Jd~ and G(p) = i (o 2 p + — 2 p _ ) for a, P 6 R. 

For investigation risk premium of the uncertainty model we give the representation of the process b t 
which is N(\fJt,pt], {0}) distributed as follows 

Lemma 2.1 For each fix t G R + , we assume that /j t be identically distributed with b\, i.e., /j, be 
N( \[i,p\ , {0}) distributed, and for t > we assume that \i t is independent from (jj ti 7J u (2 , • • • ,p. tn ) for 
each n E N and < t\, • • • ,t n <t. Let %f = ,tf , • • • ,t$} be a sequence of partitions of [Q,t], we 
define ^Q^tdt — Y.k=Q \ tnen foUtdt is N (\jJt ,pt] , {0}) distributed in (i2,L, p (i2),£). 

Proof. For t > 0, we denote a t = f {) ji t dt. For t ,s > it is easy to check that a l+s — a t is identically 
distributed with a s and independent from (a tl , a tl , ■ ■ ■ ,a tn ) for each n E N and < t\ , • • • ,t„ <t, and 
there holds 

E[a t }=jjt, —E[—a t ]=/jt. 
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Since /j, is independent from (/j f , ,/j t2 , • ■ ■ ,rt„) f° r eacn n G N and < ?i, •••,?„<?, we have 
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£[« 2 ] = £[(£^f + i-'f)) 2 ] 

jfe=0 

= £[£[(y+W-l(^-^-i)) 2 Lj*-a ft((J r ,_#)]] 

= «[(£W +1 -<f)+^-4-i)) 2 ] 

fe=0 

= Ai 2 ? 2 . (2.5) 

Forcp(x) (zCt,.iib{R) set m(? = it[<p(jt+fif)], we have 

|w(f,x) — < C|jc — y|, Vjc,yG/?, (2.6) 

where C only depend on the Lipschitz constant, and in the rest of the proof C only depend on the 
Lipschitz constant, fj and /7. 

For 8 € [0,f], since a t — 05 is independent from 05, we have 

u(t,x) = E[(p(x + a§+a t — a§)] 

= ■E'[£[<P(} ; + Of - 05) |y=.v+a 6 ]] 

= £[n(r-8,jc + a 5 )], (2.7) 

then from ( |2.5| l and ( |2.6| ), we have 

|«(f ,x) — «(f — 8,x)| 

< E[\u(t — 8,x+a§) — u(t — 8,x)\] 

< CE[\a & \]<C8. (2.8) 
Thus we prove that u(t,x) is Lipschitz in x and t . Consider the following PDE 

d,u-g(d x u) = 0, (t,x)eR + xR, (2.9) 

u\t=Q = <P, 

where 

g(oc) =Jia + —paT, Voce/?. 

To prove u(t,x) is the viscosity solution of ( |2.9| , for fixed (f,x) G /? + x we set v e C^ 2 {R + x /?) 
such that v > u and v(?,x) = u(t,x). Notice that E[a$] = /78, — £[— a§] = /j8, £[a|] = /7 2 8 2 and (2.7 ), 
by Taylor's expansion we have 

< £\y(t — 8,x+a$) — v(t,x)] 

= E [-d t v(t,x)8 + d x v(t,x)a 5 + C(8 2 + aj)] 

< -d,v(t,x)d + E[d x v(t,x)a & }+E[C(S 2 +al)} 

< -d t v(t,x)8+g(d x v(t,x))8 + C8 2 , 

from which we have d t v(t,x) — g(d x v(t,x)) < 0, we prove that u(t,x) is a viscosity subsolution of 



(2.9 1. We can prove that u( t ,x) is also a viscosity supersolution of ( |2.9| in a similar way, thus u{t,x) 
is the viscosity solution of (2.9 ). It follows thata f = J fJtdt is N(\fit,~pt],{0}) distributed. The proof 
is complete. □ 

Since the quadratic variation process < B >, of the G-Brownian motion B t is 7Y([o 2 f,0 2 f], {0}) 
distributed, from Lemma [2~T| we have 
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Lemma 2.2 <B> t is identically distributed with [q o t dt, where for each fixed t € R + O t is N( [o 2 , 2 ] , {0}) 
distributed, and for f >0o, is independent from (a tl , G t2 , ■ ■ ■ , O tn ) for each n£N and < t\ , ■ ■ ■ ,t n < 
t. 

From Lemma 1 we have that b t is identically distributed with ^fi s ds, then the price process of the 
stock can be rewritten as follows 

dS, = S, (rdt + Q,dt + dB t ), (2. 10) 

where Q, = p t — r is the risk premium which is the difference between the return rates of the stock 
and bond. It is easy to check that 8 r is N(\p. — r,p— r], {0}) distributed. 

Consider an investor with wealth Y, in the market, who can decide his invest portfolio and con- 
sumption at any time t e [0,T]. We denote % t as the amount of the wealth Y t to invest in the stock 
at time t, and C(t +h) — C(t) > as the amount of money to withdraw for consumption during 
the interval (t,t +h],h > 0. We introduce the cumulative amount of consumption Q as RCLL with 
C(0) = 0. We assume that all his decisions can only be based on the current path information £2 ( . 

Definition 2.1 A self-financing superstrategy (resp. substrategy) is a vector process (Y,n,C) (resp. 
(—Y,n,C)), where Y is the wealth process, n is the portfolio process, and C is the cumulative con- 
sumption process, such that 

dY, = rY t dt + n,dB, + %,Q,dt - dC, (2. 1 1) 

(resp. - dY, = -rY,dt + %,dB, + % t Q t dt - dQ ) (2. 12) 

where C is an increasing, right-continuous process with Co = 0. The superstrategy ( resp. substrat- 
egy) is called feasible if the constraint of nonnegative wealth holds 

Y t >0, te[Q,T}. 

3 Bid-ask pricing European contingent claim under uncertainty 

From now on we consider a European contingent claim ^ written on the stock with maturity T, here 
^ € Lq(Q.t) is nonnegative. We give definitions of superhedging (resp. subhedging) strategy and 
ask (resp. bid) price of the claim ^. 

Definition 3.1 (1) A superhedging (resp. subhedging) strategy against the European contingent 
claim t, is a feasible self-financing superstrategy (F,7t,C) (resp. substrategy (—Y,n,C)) such that 
Yt = ^ (resp. —Yt = —£,)■ We denote by ?{(£,) (resp. H' (— ^)) the class of superhedging (resp. sub- 
hedging) strategies against ^, and if 9{(Q (resp. ?('(—£,)) is nonempty, ^ is called superhedgeable 
(resp. subhedgeable). 

(2) The ask-price X(t) at time t of the superhedgeable claim £ is defined as 

X(t) = inf{x > : 3(Y,,n t ,Q) e H (£) such that Y, = x}, 

and bid-price X'(t) at time t of the subhedgeable claim t, is defined as 

X'(t) = sup{x > : 3(-Y,,n t ,C t ) e such that - Y t = -x}. 
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Under uncertainty, the market is incomplete and the superhedging (resp. subhedging) strategy of 
the claim is not unique. The definition of the ask-price X(t) implies that the ask-price X(t) is the 
minimum amount of risk for the buyer to superhedging the claim, then it is coherent measure of risk 
of all superstrategies against the claim for the buyer. The coherent risk measure of all superstrategies 
against the claim can be regard as the sublinear expectation of the claim, we have the following 
representation of bid-ask price of the claim. 

Theorem 3.1 Let ^ € Lq(Q.t) be a nonnegative European contingent claim. There exists a super- 
hedging (resp. subhedging) strategy (X,Tl,C) € ^(^) (resp. (— X' ,%,C) € !H'(— t,)) against ^ such 
thatX t (resp. X' t ) is the ask (resp. bid) price of the claim at time t. 
Let (H' s : s >t) be the deflator started at time t and satisfy 

dHl = -Hi [rds + — dB s ] , Hi = 1. (3.13) 

Then the ask-price against ^ at time t is 

X t = E[H t T Z,\£l t ], 

and the bid-price against ^ at time t is 

X' t =-E[-H' T ^], 
Proof. By G-Ito's formula we can check that 



H,=exp{-[ rds+ —dB s + - ( — yd<B>,]} (3.14) 

JO CTj 2JoO s 



is the solution of (3.13 ). Define the stochastic process X from 

H t X t =E[H T QCl t }=M„ 

then M t is a G martingale. By the G-martingale representation theorem (|25|), there exists unique 
decompasition of H,X t as follows 

H t X t =E{H T £,}+ f'$ s dB s -K t , 
Jo 

where {p r } e Hq(0,T), {K t } is a continuous, increasing process with Kq — 0, and {— K t }o<t<r is 
a G-martingale. Set n t = [H,~ 1 p, + X, |-] . Then H,X, = E [H T $] + Jo H s (jt s - X s )dB s - K, . Define 
C(t) = f H~ K s ds, then C(f) is nonnegative and increasing process with C(0) = 0. We prove that 
(E[H' T Z,\Q. t ],K t ,C t ) G is a superhedging strategy against t,. 

For any superhedging strategy (Y t ,%,Q) against ^, by G-Ito's formula and a t dt = d < B > t 
which given by Lemma 2.2 we have 



H t Y t =H T ^- f H s {% s -Y s — )dB s + f H s dC s , 

Jt O s Jt 



Taking the condition G-expectation on both side of (3.15 I with respect to Q. t , notice that C t is a 
nonnegative right continuous process and E[J t T H s (n s — Y s Q s )dB s \£l t ] = 0, we have 

H t Y{t)>E\H T ^\0. t ], 

i.e. 

Y(t)>E[H , T ^\Q. t ]=X„ 

which prove that X t = E[Hjt,\Q.,] is the ask price against the claim t, at time t. Similarly we can 
prove that X\ = —E[—H t T £,\Q. t ] is the bid price against the claim £ at time t . □ 
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4 G martingale transform and bid-ask dynamic pricing mecha- 
nisms 

In this section, we will construct a time consisten G expectation E G , and transfer the uncertainty 
model into the sublinear space (£l,L G (Q.),E G ) which is correspond with a sequence of the risk- 
neutral probability measure space. 

Define a sublinear function G( ) as follows 

G(oc) = i(a 2 oc + -g 2 or), Voce/?. (4.15) 

For given <p £ Cb.u p (R), we denote u(t,x) as the viscosity solution of the following G-heat equation 

d t u-G(d xx u) = 0, {t,x) £(0,oo)xR (4.16) 
u(0,x) = (p(x). 



Remark 2 The G-heat equation (4.16 1 is a special kind of Hamilton- J acobi-Belhnan eqaution, also 
the Barenblatt equation except the case O = (see and The existence and uniqueness 

of (4.161 in the sense of viscosity solution can be found in, for example M3V , MOV , and M7V for 
C l,i -solution if O > 0. 

Denote 



B,=B,-rt, (4.17) 

the stochastic path information of {B t }t>o up to t is the same as {B t },>Q, without loss of generality 
we still denote £l t as the path information of {B t }t>o up to t. For CO £ Q. consider the process 
(B t )((o)=(O t ,t£ [0,°o), we define E G [] :L ip (£l) — >Ras 

E G [( ? (B t )]:=u(t,0), 

and for each s , t > and t\ , ■ ■ ■ , £ [0, t ] 

E G [i?(B tl , ■ ■ ■ ,B lN ,S t+s -B t )\ := E G [\]f(B tl ,S, N )} 

where ,•• • ,x N ) = E G {<p(x u - ■ ■ ,x N ,B s )]. 

For < ti < ti < ■ ■ ■ < t{ < ti+i < • • • < tu < +°°, we define condition G expectation with respect 
to £l tj as 

E G [<?(Bt l ,Bt 2 ,B tj+l -B tj ,- - ,B, N — B tN l ) | £1^ ] 

:= y(B tl ,B t2 —B tl ,- - ,B t . — B tj _j), 

where ,■• • ,*,-) = E G [q>(xi , • • • ,x u B tM —B tj ,- - ,B, N ~B, N _ l ). 

We consistently define a sublinear expectation E G on L,„(f2). Under the sublinear expecta- 
tion E G we define above, the corresponding process (B t )t>o is a G-Brownian motion and B, is 
A^({0}, [g 2 f ,CT 2 ?]) distributed. We call E G as G-expectation on Lf P (£2). 

Denote L p G (Q.),p > 1 as the completion of L,- p (£2) under the norm ||X|| P = (£' G [|X| / ']) 1 / /J , and 
similarly we can define L p G (£l t ). E G can be continuously extended to the space (Q.,L G (£l)), From 
now on we will work in the sublinear expectation space (Q.,L G (Q.),E G ). 
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The price dynamic process of the stock (2.4) can be rewritten as follows 



dS t =S,(rdt + dB t ). (4.18) 

Denote D(t) := e~ rt be the discounted factor, with the discounted processes Y t = D(t)Y t , % = 
D(t)n t , and dC t — D(t)dC t , using G-Ito's formula, we can write the self-financing superstrategy 
(Y,n,C) (resp. substrategy (— Y,K,C)) satisfying 

df, = H t dB, - dC. 

(resp. — dY, = K t dB t — dC ). 
The superhedging (resp. subhedging) strategies and ask (resp. bid) price of the claim ^ which we 



defined in Definition |3.1| can also be characterized using discounted quantities. 

For the nonnegative European contingent claim "t, € L g (Q.t), we define G dynamic pricing mech- 
anism for the claim t, as follows 

Definition 4.1 For t £E [0, T\, we define G dynamic pricing mechanism as Ef T : L G (£lj) — >■ L G (Q. t ) 

E^ T [-}=E G m t }. 



By the comparison theorem of the G-heat equation (4.16 1 and the sublinear property of the 
function G(-), the G dynamic pricing mechanisms have the following properties 

Proposition 4.1 For t e [0, T] and ^i,^e L 2 G (Q. T ) 

(i) E? T [^]>E? T M,if^>Z, 2 ; 

(ii) EG T g]=Z > ; 

(iii) E^+^ 2 ]<E^U+E^ T U; 

(iv) Ef T [A£] = lEf T g] for X > 0; 

(v) E° [Efr &}] = Ef T g] for0<s<t. 

Theorem 4.1 Assume that ^ = (])(S:r) € L G (£lj) be a nonnegative Eur opean contingent claim, and 
E^ T [•] be the G pricing dynamic mechanisms defined in Definition 4. 1 The ask price and bid price 
against the contingent claim ^ at time t are 

M a (f,S ( )=e- r ( r -')£, G r g] and u h (t,S t ) = -e- riT -'^E^ T [-^} (4.19) 

respectively. 

Proof. It is easy to check that {uf }o<;<r satisfying 

D(tX=E^ T [D(Tn 

then M t = D(t)uf is a G-martingale. By a similar way we used in the proof of Theorem 3.1 we can 
complete the proof. □ 

Remark 3 In the Properties ^. 1| (Hi) and (iv) imply that 

+ (1 - a)fc] < oE^gi] + (1 - a)E^ T g 2 ]Jor a e [0, 1], 

which means that the G dynamic pricing mechanism is a convex pricing mechanism, (v) means the 
G dynamic pricing mechanism is a time consistent pricing mechanism, under the G dynamic pricing 
mechanism the ask price u a (s,S s ) = e~ r ( T ~ s } E^ T g] (resp. the bid price u b (s,S s ) = — e^ r ^ T ^Ef T [— Q) 
at time s (s <t <T) against the claim ^ with maturity T could be regards as the ask ( resp. bid) price 
at time s against the claim u"(t,S t ) — e~ r ( T ~^Ef* T [^] (resp. u b (t,S t ) = — e~ r ' T ~ t ^ E^ T [— £,]) with 
maturity t. 
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5 Markovian case 



We assume that the stock price dynamic satisfying the following SDE (f > 0) 

dS'/ x = S' s x {rdt + dB t ), se[t,T], (5.20) 
Si"" = x. 



4.1 



its 



For given a nonnegative European contingent claim £ = §(St) G Lq(£2), from Theorem 
ask price and bid price at time t are 

u b (t,x) = -e- r ^E^ T [-^)], 

we establish the relation between ask (resp. bid) price with the viscosity solution of a full nonlinear 
PDE. 



Theorem 5.1 Assume that the stock price dynamic process satisfying (5.20), \ = §{S't) G Lq(CIt) 
be a nonnegative European contingent claim written on the stock with maturity T, and (j) : R — > R 
be a given Lipschitz function. The ask price of the contingent claim u a (t,x) = e^^^'E^ [^] is the 
viscosity solution of the following nonlinear PDE 

d t u a (t,x) + rxd x u a {t,x) + G(x 2 d xx u a (t,x))-ru a (t,x) = 0, (f ,x) € [0, T) x R (5.21) 

u a (T,x) = iSf{x). 

The bid price of the contingent claim u b (t,x) = —e^ r ^ T ^Ef T [— is the viscosity solution of the 
following nonlinear PDE 

d t u b (t,x)+rxd x u b (t,x)-G(-x 2 d xx u b (t,x))-ru b (t,x) = 0, (t ,x) G [0, T) x R 

u h (T,x) = ty(x). (5.22) 

Proof. With the assumption that "% G Lq(Q.t), Peng in G3l prove that v(t,x) =E^ T [^(S'j)] satisfying 
|v(*,Jt) - v(?,*')| < C\x-x'\, \v{t,x)\ < C(l + |*|), 

and for 8 € [0, T - 1] 

|v(f,jc)-v(* + 6,jt)| <C(l + |x|)(5'/ 2 + 8), 5e[0,T-t], 

v(t,x)=E G [v(t + 8,!%J), 

where C is only dependent on the Lipschitz constant. 
We can easily get that 

\u"(t,x)-u a (t,x')\ <C\x-x'\, \u a (t,x)\ <C(1 + |*|), 
\u a {t,x) - u a (t + 8,*)| < C(l + |x|)(8 1/2 + 8), 

and 

K a (M)=£ G [e~ r V(* + S,s£ 5 )]. 
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For fixed (t,x) € (0, T) x R, let \|/ € C^ ,3 ([0,r] x fl) be such that y > u" and = u"(t,x). 

By Taylor's expansion, we have for 8 G (0, T — t ) 

< £ G [ e - r8 ¥ (f + o,^ 5 )- ¥ (f,*)] 



< -£ a [x z a„ ¥ (f ! x)(<B> ;+8 -<B> r )] 

+(d t y(t,x) + rxd x y(t,x) - ry(t,x))8+C(l + \x\ + |x| 2 

< (d,\\l(t ,x) + rxd^ (t ,x) + G(x 2 d xx \\l(t ,x)) - r\\t(t ,x))d 



! )8 3 



-C(l 



! )5 3/2 , 



for 8 4- 0, we have 



3(\(/(f,x)+rx3, : \|/(f,x) + G(jc 2 3.«\|/(f,x))-r\|/(f,x) > 0, 



which implies that u"(t.x) is the subsolution of the nonlinear PDE (5.21 ), and by a similar way 
we can prove that u a (t,x) is the supersolution of ( |5.21|). Thus, we prove that the ask price against 
the claim ^ at time t is the viscosity solution of ( 15.211) . Similarly, we can prove that u h (t,x) = 
-e- r{T -^E^ T [^{S'j X )] is the viscosity solution of (15.221). □ 
Assume that the stoke price solve the following STJET 



d\ogS' u 



rdt + d t dW tl u e [t,T] 



(5.23) 



where (W r ) r >o be a 1-dimensional standard Brownian motion defined on a probability space (£2, J- ,P) 
and the filtration generated by W t is {% '■ < t < T). Assume that (c t )t>o is an adapted process such 
that a t E \g,g]. It is well known that the Black-Scholes type price against the European contingent 
claim i; = (\>(S f j) at time t is u(t,x) = e~ r ( T ~ t " > E[ty(S t j x )\y r t ] which is the solution of the following 
PDE 



d t u(t,x) + rxd x u(t,x) + -x 2 afd xx u(t,x) — ru(t,x) 

u(T,x) 



0, (t lX )eR + xR (5.24) 

4»(x). 



Corollary 5.1 Assume that § : R — > R be a given Lipschitz function, u(t,x) is the Black-Scholes type 
price against the claim ^(Sj X ) which satisfies llblack). Assume that u b (t,x) and u a (t,x) satisfying 
(5.22)) and (5.21 ), i.e., be the bid price and ask price against the European contingent claim ^ = 



§{Sf ) at time t. Then we have 



u b (t,x) < u(t,x) < u"(t,x). 



(5.25) 



Proof. By using the comparison theorem proposed by Peng in 
□ 



we can easily prove the Lemma. 



e 



Lemma 5.1 Assume that the price process (S s ' x ) s > t of the stock satisfies (5.20), \ = 
L 2 j(Q.f) be a nonnegative European contingent claim written on stock with maturity T, and (j) : 
R — > R be a given Lipschitz function. 

(I) If (])(•) is convex (resp. concave ), for any t G [0,T] the ask price function u"(t,-) is convex 
(resp. is concave), and u a (t,x) satisfies (5.21 ) with G(x 2 d xx u") — ^a 2 x 2 d xx u a (resp. — ^a 2 x 2 d xx u a ). 
If (()(•) is convex (resp. concave) on interval (a,b) € R (a < b) the ask price function u a (t,-) is 
convex (resp. concave) on (a,b) for any t G [0, T], 
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(II) If §(■) is convex (resp. concave), for any t G [0,7*] the bid price function u b (t,-) is con- 
vex (resp. is concave ), and u b (t,x) satisfies (5.22) with G(—x 2 d xx u b ) = — ^0 2 x 2 d xx u b (resp. = 
~^0 2 x 2 d xx u b ). Ifty(-) is convex (resp. concave) on interval (a,b) G R (a < b) the bid price function 
u b (t, •) is convex (resp. concave) on (a,b)for any t 6 [0, T\. 

Proof. We only prove (I), the proof of (II) is similar with the proof of (I). 
1. First we prove that u a (t,x) is convex if §(x) is convex. 

For x\,X2 G R and for given a G [0, 1], by G-Ito's formula, we get the stock price process as 
follows 

c f,ou"i+(l-a)x 2 

= (axi + (1 - a)x 2 ) exp (r(T - 1) + B T - B, - ~ (< B > T - < B >,)) 
= <xxie\p{r{T -t)+B T -B t -^(<B> T - <fi >,)) 

+(1 - a)x 2 exp (r{T -t)+B T -B t - l -(<B> T -<B >,)) 



aS/ 1 



?/ 2 . 



(1 -a)Sj 

Since the G pricing dynamic mechanism is a convex pricing mechanism, we have that 

u"(t, oui + (1 - ate) 
= e-^E^iS'^ 1 -^ 2 )] 
= e- r[T -^Ef T [§(aS' T X[ + ( 1 - a)S'j X1 )] 



(5.26) 



< 



-r(T-t) F G 



a^(5'/') + (l-aH(5 f / 2 )] 



< ae-^E^ T feCS*/ 1 )] + ( 1 - a)e-^ T ^E^ T )] 

= Uu a (t,X l ) + (l-u)u"(t,X2), 



(5.27) 



which prove that u a (t,x) is convex on R, d xx u a is nonnegative on R and G(x 2 d xx u a ) = ^a 2 x 2 d xx u". 

If (|)(x) is convex in interval (a,b), for x\ ,x% G (a,b), (5.26 ), (5.27 ) hold on (a,b) and u"(t ,x) is 
convex on (a,b). 

2. We will prove that u"(t,x) is concave if (j)(x) is concave. 

For given process (c,),>o such that a, G [o,a], we denote «£( f tne solution of the following 
PDE 



1 



dtU a a (t,x) + rxd x u%(t,x) + -O t x d xx u^(t,x) 

u«(T,x) 



0, (t 1 x)eR + xR 



then u a a (t,x) =E t [$(xexp[r(T -t)+o,(W T - W t )- \o, (T-t)])], where W t is the standard Brownian 
motion defined on a probability space (£1, f , (!ft),P), and E t is the corresponding condition expec- 
tation. Denote S/ a = xexp[r(T — f) + o t {Wr — W,) — jO?(r — ?)]> for jci G R if is concave 
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on R, we have 

u%(t,ttx\ + (1 -a)x2) 
= E t [if{{ax x + (1 - a)x 2 )exp[r(r -t) +O t (W T - W t ) - ~o}(T -?)])] 

= £^(a^ + (l-a)S^)] 
> E t [a^S'^) + (l-aMS'^)} 
= aE t MS'^)} + (l-a)EMS^)} 

= au%(t 1 xi) + (\-a)u a a {t,x 2 ), (5.28) 
which mean u%(t, •) is concave on R, and the function U (t,x) = sup OG r og i u a a {t ,x) is concave on R, 



U{t,ax\ + (1 -a)x?) > at/(*,xi) + (1 -a)i/(f,*2)- 



(5.29) 



Since the operator sup OG [ aS ]£' f is a convex operator and U(t,x) — sup ae ^ a ^E,[<^(Sj X a )], using the 
similar argument in Theorem 5.1 we can prove that U(t,x) is the viscosity solution of the following 
HJB equation 

d,U(t,x) + rxd x U(t,x) + l > sup [a 2 x 2 d xx U(t ,x)] = 0, (t,x) e R + xR 



U{T,x) 



which can be rewritten as 



3 f f/(f ! x) + rx3 x f/(f,x)+G(x 2 a xx f/(f,x)) = 0, (t,x)eR + xR 

U(T,x) = <|>(jt). 



(5.30) 



In |23i, Peng prove that U{t,x) — E^ T [<^(Sj X )] is the unique viscosity solution of (5.30 1, from (5.29 1 
the solution U (t,x) is concave for x on R. Thus we prove that u a (t ,x) — e~ r ( T ~')U(t,x) is concave 
forx on R and G(x 2 d xx u") = la 2 x 2 d xx ii" . 



If (j)(x) is concave on (a,b), ( |5.28 I and (5.29 1 hold on (a,b), which prove that u"(t,x) is concave 
on (a,b). We finish the proof of (I). □ 



6 Monotone characteristic finite difference schemes 

In this section we will consider numerical schemes for the nonlinear PDE ( |5.21[ ) (similar for ( |5.22| >). 

6.1 Characteristic finite difference schemes 



Define u(z,x) = u"(T —t,x), the nonlinear PDE (5.21 1 (similar for (5.22 1) can be written as 

d x u(x,x) — rxd x u(x,x) — x 2 G(d xx ii(i,x)) + ru(x,x) = 0, (t,x) 6 [0, T) x R, (6.1) 

u(0,x) = §(x). 



First, we consider the boundary conditions of (6.1 1. As x —> 0, equation (6.1 1 becomes 

3t«|a-=o = -ru(x,0). 



(6.2) 
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For x — > oo, normally for sufficient big enough x there holds d xx u ~ 0. We set x = S max with 
Smax big enough such that the payoff of the contingent claim has asymptotic form. We consider the 
Dirichlet condition as follows 

u\x=S am =g(l,Smi*), (6.3) 

where g(.,l) can be determined by financial reasoning for some given contingent claim, normally in 
the following asymptotic form 

g(T,S m ax) = b(x)S max +c(x). 

We assume that b(x),c(x) are bounded such that 

|wMmax)|<Q,, (6.4) 

where Q, is a constant. 



There is convection term d z u{x,x) — rxd x ii(x,x) in (6.1 ), if the convection dominate the diffu 



sion finite difference discretization for (6.1 ) could leads numerical oscillation, we consider discrete 
the convection term along the characteristic direction (|9|). Denote w(x) = [1 + r 2 x 2 ] l l 2 , the direc- 
tion derivative along characteristic direction is ^ = — rx §^)- (6.1 )-( 6.3 ) is equivalent to the 
following PDE 

¥ (x)|-su P$? < o < g i: o M = 0, (z,x) e (o,t] x (o,s raax ), 

3x«|x=o = -ru(x,0), u\ x=Sma = g(l,S mia ), (6.5) 
u(0,x) = cb (jc) , 

where LPu — ^x 2 d xx u — ru. 

Now we define space partition of / = [0, S max ] . Let / = [0, S max ] be divided into N sub-intervals 

/, = (xi,Xi+i), i = 0,...,N-l, 

with = xq < x\ < ■■■ < Xff = S max . For each i = 0, ...,N — 1, let Ax,- — X{+\ — Xj, Let {ti}^L be a 
set of partition point in [0, T] satisfying = to < t\ < ■ ■ ■ < tm = T, and denote Af„ =t n — t n -\ > 0, 
where M > 1 is a positive integer. 

Let uf be a discrete approximation to u(t n ,xi) . Denote I" = x, + rx[At n+ \ , for At n+ \ small enough 
such that x" € Denote u" be a discrete approximation to u(t n ,x"), here we define u" as the 



linear interpolate function of u". The implicit characteristic finite difference scheme for (6.1 1 is as 
follows: 



For boundary x = 



it 



n+1 n 

_ „„n+l 



for i = 1,2,- •• ,N-\ 



rul +i , (6.6) 



„«+l ,.n , 

"t= sup [(L / r + " n+1 ) I -]] + ^, (6-7) 

where (L„ u n )i denotes the central difference discrete of {LPu) at note (t„,Xj), i.e. 

(Lf u ")i = a?WK_i + P?(c?K +1 - (a5»(o?) + p?(o?) + rK, (6.8) 
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where a" and p" are defined as follows 

(o?) 2 xf 



a?(o?)=/ 

[Xi—Xi-l)[Xi+l —Xi-l) 



*\ 2v 2 (6-9) 



(jCj+l -Xi)(Xi+l -Xi-l) 

It is easy to check that a" , p" > 0. We define 

U" = [u'q , u" , • • • , u" N _ j , Uff] T , 
fl" — \ii n u" ... if u n ] T 

a n = [o1,-,o n N _ l ] T , 

and 

(A"(a")U")i = o?MWi-x +P?(o?)«? +1 - (a?(o?) + p?(of) + rK. 
For notational consistency, we denote SJ] = l+r \j — -Uq, — m^ +1 , and enforce the first row and the 



last row of A to be zero. Then the discretization scheme ( 6.7 1 can be write as the following equivalent 
matrix form 



[I-At n+l A"+ l ]U"+ l =U" 

ar 1 =argsu Po „ +le{ ^ } {(A»+ 1 f/"+ 1 ),} (6 - 10) 



where A n+l = A(6 n+1 ) andd" +1 =c"+ l (U" +l ). 

It is easy to see that — A" +1 has nonpositive off-diagonals, positive diagonal, and is diagonally 
dominant. We have the following theorem 

Theorem 6.1 Matrices -A" +1 and I - Af" +1 A" +1 are M-matrices ([26]). 



The equation ( |6.1| l has unique viscosity solution, and satisfies the strong comparison property 
(see J6) and H0|), then a numerical scheme converges to the viscosity solution if the method is 
consistent, stable and monotone. 

Let h = max {Ax, A?} be the mesh parameter, where Ax = max,Ax/,Af = max„ Af„. Assume that 
the partition is quasi-uniform, i.e., 3Ci,C% > independent of h, such that 



Cih< Axi,At n <C 2 h (6.11) 



for < i < N- 1 and 1 < n < M. 
We denote 



G? + \h,ul^,u? +l ,u£l,i%) = sup (6.12) 

Ar»+1 a"+'e{o.o} Af »+1 

then the discrete equation at each node can be written as the following form 

G;' +1 (/*,<_VX + VmX)=0 (6.13) 

Lemma 6.1 (Stability) The discretizaition ( |6.10p is stable i.e. 

<max(||C7°||.,C fr ). (6.14) 
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Proof. The discrete equations are 



(1 +Af„ +1 (a; ,+1 +p? +1 +r)K +1 = u n i +Af„ +1 < +1 <+ 1 +Af„ +1 pf '<+/. 
Since we take u" as linear interpolation of u" and m" +1 in the discretization, there holds 

(l+Ar„ +1 « +1 + p;. ,+1 +r))|< +1 | < ||f/"|| 00 + Ar„ +1 « +1 +p^ 1 )||C/" +1 | 



If lie/ 



n+l|| —\„n+\\ 



< j <N, then we have 



(l+At„ +1 (af 1 + p? +1 + r))||C/"+ 1 ||« J <||C7' ! ||« J +A?„ + i(af 1 + pf 1 )||C/" +1 | 
which implies that 

If j = or j = N, then ||oo = |"o +1 1 < \u n \ or \\U" +1 \\„ = \u" N +1 \ < C b . 



Thus, we have 



which complete the prove. □ 



\U 



«+ii 



<max(||t/ || ro ,C fo ) 



Lemma 6.2 ( Consistency) For any smooth function v with v" = v(f",x,'), the discrete scheme (6. 10 1 
is consistent. 

Proof. Using Taylor's expansion, we have 

|(iaa„«-rv)7-(Ljv"),|=0(Ax), 



and using expansion along characteristic direction c 



,,«+! 



y dc H At, 



n+l 



= 0(Ax + At). 



where v" = v(t n ,x"). 

For smooth v, by Taylor's expansion, we can derive the discretization error as follows 



.3v. 



0P^)? +1 - sup (L 



a,, n+l \ 



8x 



og{o,o} 

v" +1 -v", 



v" +1 -v? 



Af, 



-i- sup (Z° v" +1 ), 
n+l a" +1 6{a,5} 



Ar 



n+l oe{a,o} 



sup |(LV' +1 ),-(L°v 



/i 1 



= 0(At+Ax) 
which prove the consistency of the discretizaiton schem. □ 
Lemma 6.3 (Monotonicity) The discretization ( |6.10| ) is monotone. 
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Proof. For ; = or i = N the Lemma is trivially true. For < i < N, we write equation (6.12 1 in 
component form 



Af„ + i Af„ + i 

- sup [a; ,+1 (a? +1 x+l + p; ,+1 (c? +1 )<+ 1 1 - (a; ,+1 (a^ 1 ) + ) + r)< +1 ] . 

a" +1 e[a.a] 

For e > 0, we have 



; sup {-p; ,+1 (o; ,+1 )e} = -e inf _{p; ,+1 (o" +1 )}<0. 

a"+'G[a,S] a"+'G[o : o] 
With the similar argument we derive 



Cr, (tl, U i x + e, Uj , , M ( - ) - Lr i (ll, U i x , U i , U i+1 , U { j < U. 



It is easy to check that 



r" 1 + v lh ,, n + l ,, n + l n n j-c\ f n +^lu r/"+l r;""! 



<0. 



+ 1 



Thus we prove the discrete scheme (6.IO1 is monotone. □ 



The discrete scheme (6.IO1 is consistent, stable and monotone, from [5|, we have the following 
convergence theorem 



Theorem 6.2 (Convergence to the viscosity solution) The solution of the discrete scheme (6.101 
converges to the viscosity solution of equation (|6.1|). 



6.2 Iterativ Solution of Discrete Algebraic System 



In the previous subsection, we have show that the solution of the discretization (6.7 1 converges to the 



viscosity solution of the nonlinear PDE (6.1 1. Since the implicit scheme leads a nonlinear algebraic 



system (6.IO1 at each timestep, the discretization is not a practical scheme. In this section, we aim 
to solve this discrete scheme by a practical iterative method. 



Iterative Algorithm for ( 6. 10 1 

1. n = 

2. Set k = and u k = U" 

3. For& = 0,1,2,- •• 
Solve 



[I-At„ +l A" +l (a k )}u k+l =0" 
ofeargsup^^jl^+Ha)^],} 



(6.15) 



4. Ifmax,- '"' u ! +l 

mz\x{scale,Uj 

5. Set U" +1 



< tolerance then quit, else k = k + 1 go to 3. 



u k +\a n+1 



6 and n 



1, go to 2. 



The term scale is used to ensure that unrealistic levels of accuracy are not required when the 
value is very small. In the iterative algorithm a is given by 



a, if 
g, if 



(u k +l ~U k )/{x i+l -^Q-^f-HjLJ/^.-.V,-! ) 

{u k + , - u k ) / (x i+ !-*,■)-(«*-«*_,)/ (Xj 1 ) 
X i+ l-Xi-l 



>o 

<0 
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Theorem 6.3 (Convergence of the Iterative Algo rithm) The iteration algorithm (6.15 I for (6.101 
converges to the unique solution of equation (6.7 \ for any initial iterate u°. 

Proof. First, we will prove that ||m ||«, is bounded independent of iteration k with a similar argument 
we use in Lemma |6~T| 

We can write ( |6.15[ ) in component form as follows 

( 1 + Ar„ +1 «+> (6*) + p; i+1 (5*) + r))u k+1 - Af„ +1 < +1 (a*)^ 1 - Ar„ +1 pf 1 {&)u\ %\ = u'} . 



Then 



1,^+1 1 



< 
< 



(1 + AWi (a; ,+1 (3*) + p; ,+1 (6*) + r))\uf 

„n+l (zk\ -k+\ i a, - (o*)^ 4 ■ 

sl n+\{^k\ i dm+1 I z-k\\\\ ~k+\ I 



|«"| + |Af„ +1 af +1 (^)«r/ +AWiPr (^K+i 



Thus 



If/" I 



-Ai B+ i(a? +1 (a*) + p? +1 (e*))||fi* 



<||C7"|U<max(||C/°|U,a), 



which means that \\u + \\ m is bounded independent of k. 

Now we will prove that the iterates {u k } form a nondecreasing sequence. From (6.15 1, the 
iterates difference u k+l — u satisfy 

[/-Af„ +1 A" +1 (d* : )](M i:+1 -fi A: )=Af„ +1 [A" +1 (d' : )-A" +1 (d' : - 1 )]M A: . 



(6.16) 



Notice that 



of G argmax OGfeg} {[A" +1 (°)]«1,-}, 



1,1*1 



the right hand side of (6.16 1 is nonnegative, i.e. 

At n +i{[A" +l (a k ) -A n+l {o^ 1 )}^}, >0. 

Consequently, 

[{I-te n+l A n+ \&)){u k+l -u k )]i > 0. 
we know that matrix / — Af„ + iA" +1 (a*) is an M-matrix and hence 
[I-At^A^ia'T^O. 



From Theorem 



6.1 



(6.17) 
(6.18) 



From ( |6.17| i ( |6.18| l, we can derive that 



u k+l -u k >0 1 



which prove that the iterates form a nondecreasing sequence. The iterates sequence {u } is nonde 



creasing and bounded, thus the sequence converges to a solution, i.e., 3u such that \\u k — u\\ 
and 

[(/-Ar„ +1 A" +1 (d)]fi = S" 
en e argsup aGtes} {[A" +1 (a)w] ! } 



oo^O 

(6.19) 



Since the matrix / — Af„ + iA" +1 (d) is a M-matrix, thus the solution of (6.19 1 is unique which finish 
the proof. □ 

In next section we will simulate the ask (resp. bid) price of the contingent claim by using mono- 
tone characteristic finite difference schemes for (5.21 1 (resp. (|5.22|i). The convergence of the solu- 



tion of the monotone characteristic finite difference schemes (|6.7|i to the viscosity solution of (6.1 



guarantee the simulation ask (resp. bid) price of the contingent claim convergence to the correct 
financial relevant solution. 
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7 Examples and simulations 



In this section, we will give simulations for the bid-ask pricing mechanisms of contingent claims 
under uncertainty with payoff given by some functi on §( St)- In computer simulations, we only 



make the numerical program for the nonlinear PDE (5.21 1, since the bid price u (t,x) := —u(t,x 



where u(t,x) is the viscosity solution of (5.21 1 with terminal condition u(t,x) = — 
Example 7.1 Digital call option under uncertain volatility. 
We consider a digital call option with the payoff as follows 



<t>(Si 



1, S T >K 
0, S T < K 



The strike price of the digital option is K = 100 and the maturity is six months T = 0.5. The volatility 
bounds are given by a = 0.15,(7 = 0.25 and short interest rate is r = 0.10. 

We use the numerical schemes constructed in Section 3 to compute the nonlinear PDE (5.21 1 
with the payoff function as initial condition and the boundary condition is (^(Smax) = 1. we choose 
the grid as As — I, At = 0.0025 and iterative tolerance= 10~ 6 . We plot Fig. [lithe ask (top left) pric 
and bid (top right) price surfaces of the digital call option. 

Example 7.2 Butterfly option under uncertain volatility. 

The second example is a butterfly option with the payoff as follows 

§(S T ) = max{S-K u 0)-2ma\(S-{K 1 +K 2 )/2,0)+max(S-K 27 0), 

the boundary condition is M'^Smax) —0,Ki = 90 and K 2 = 110. The other parameters used here 
are the same as that we used in example 1 . The ask price (down left) and the bid price (down right) 
surfaces of the butterfly option are shown in Fig. [Tj 

Fig. [Tjshows that the ask price surface is above the bid price surface for the both claims, and the 
ask (resp. bid) price dynamic keeps the monotone intervals and convex (resp. concave) intervals of 
the corresponding payoff function which verify the theoretical results we showed in Lemma |5~T| 
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